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Summary 


Swine acute diarrhea syndrome coronavirus (SADS-CoV), a novel coronavirus, was first 
discovered in southern China in January 2017 and caused a large scale outbreak of fatal 
diarrheal disease in piglets. Here, we conducted a retrospective investigation of 236 samples 
from 45 swine farms with a clinical history of diarrheal disease to evaluate the emergence 
and the distribution of SADS-CoV in pigs in China. Our results suggest that SADS-CoV has 
emerged in China at least since August 2016. Meanwhile, we detected a prevalence of 
SADS-CoV (43.53%), porcine deltacoronavirus (8.83%), porcine epidemic diarrhea virus 
(PEDV) (78.25%), rotavirus (21.77%) and transmissible gastroenteritis virus (0%), and we 
also found the co-infection of SADS-CoV and PEDV occurred most frequently with the rate 
of 17.65%. We screened and obtained two new complete genomes, five N and five S genes of 
SADS-CoV. Phylogenetic analysis based on these sequences revealed that all SADS-CoV 
sequences in this study clustered with previously reported SADS-CoV strains to formed a 
well defined branch that grouped with the bat coronavirus HKU2 strains. This study is the 
first retrospective investigation for SADS-CoV and provides the epidemiological information 
of this new virus in China, which highlights the urgency to develop effective measures to 


control SADS-CoV. 
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1 | INTRODUCTION 


Swine acute diarrhea syndrome coronavirus (SADS-CoV) is a newly discovered coronavirus 
which is an enveloped, positive and single-stranded sense RNA virus with a genome size of 
approximately 27kb (Gong et al., 2017; Pan et al., 2017; Zhou et al., 2018). SADS-CoV 
belongs to the family Coronaviridae which contains four genera, Alphacoronavirus, 
Betacoronavirus, Gammacoronavirus and Deltacoronavirus (Woo et al., 2012; Woo et al., 
2010). So far, six coronaviruses have been identified from pigs, which include porcine 
epidemic diarrhea virus (PEDV), porcine respiratory coronavirus (PRCV), SADS-CoV and 
transmissible gastroenteritis virus (TGEV) that all belong to the Alphacoronavirus genus, as 
well as one betacoronavirus, porcine hemagglutinating encephalomyelitis virus (PHEV) and 
one deltacoronavirus, porcine deltacoronavirus (PDCoV) (Lin et al., 2016; Woo et al., 2010; 
Wesley at al., 1991). Among these viruses, SADS-CoV is the most newly discovered 
coronavirus, which has been first reported in 2017 in China and is considered to be an 
HKU2-related coronavirus with a bat-origin (Gong et al., 2017; Zhou et al., 2018). In January 
2017, SADS-CoV was detected in a swine farm and subsequently spread rapidly to three 
other farms in Guangdong Province and caused the fatal swine acute diarrhea syndrome 


(SADS) characterized by the clinical signs with severe, acute diarrhea and rapid weight loss 
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of piglets. The symptoms of SADS-CoV are similar to those that caused by other swine 
enteric coronaviruses such as PDCoV and PEDV, but SADS-CoV is more harmful than these 
viruses because it has led to the death of almost 25,000 piglets in a short time and resulted in 
more significant economic losses (Sun et al., 2016; Dong et al., 2015; Zhou et al., 2018). So, 
it is urgent to investigate the molecular epidemiology and transmission patterns of 


SADS-CoV for establishing effective controls for this new coronavirus. 


In the present study, we performed the retrospective PCR testing on diarrheal samples 
from 45 swine farms in Guangdong Province to evaluate the emergence and the distribution 
of SADS-CoV in pigs in China. The prevalence and co-infection information of SADS-CoV 
from eleven SADS-CoV-positive farms was provided. The sequences of SADS-CoV, 
including two complete genomes, five nucleocapsid protein (N) genes and five spike protein 
(S) genes, were also identified and characterized to investigate the phylogenetic relationships 


of SADS-CoV. 


2| MATERIALS AND METHODS 


2.1 | Sample collection 


A total of 236 clinical diarrhea samples including feces and intestinal contents were collected 
from piglets and sows from 45 swine farms in the cities Qingyuan and Shaoguan of North 
Guangdong Province between August 2016 and May 2017 in accordance with the 


recommendations of National Standards for Laboratory Animals of the People's Republic of 
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China (GB 149258-2010). Samples were preserved at -80 °C from the time of original receipt 


until use. 


2.2 | Nucleic acid extraction and molecular diagnosis 


Samples were homogenized in phosphate-buffered saline (PBS) (20 % w/v), frozen and 
thawed three times, then centrifuged for 10 min at 10,000xg. Viral nucleic acid was extracted 
following the manufacturer's recommendations of AxyPrep™ Body Fluid Viral DNA/RNA 
Miniprep Kit (Axygen Scientific, Inc). The virus nucleic acid was stored at -80°C until PCR 
was performed. A pair of primers (forward primer 
5'-GGTCCCTGTGACCGAAGTTTTAG-3’, reverse primer 5'- 
GCGTTCTGCGATAAAGCTTAAAACTATTA-3) was designed to detected SADS-CoV 
based on the conserved N gene of this virus. One step RT-PCR using PrimeScript" One Step 
RT-PCR Kit Ver.2 with Dye Plus (Takara, Biotechnology, Dalian, China) was carried out to 
amplify the target fragments by the following thermal profile of 50°C for 30min, 94°C for 
3min, 35 cycles of denaturation at 94°C for 30 s, annealing at 55°C for 30 s, an extension at 
72°C for 30 s and a final step of 72°C for 5 min. Four other diarrheal pathogens including 
PEDV, PDCoV, rotavirus (RV), and TGEV from SADS-CoV-positive farms were also 
tested by RT-PCR according to the previously described methods (Mai et al., 2017; Liu et al., 


2015; Stevenson et al., 2013). 
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2.3 | Amplification of the N gene, the S gene and the complete genome of SADS-CoV 


Specific primer pairs based on reported SADS-CoV strains (GenBank accession numbers: 
MF094681-MF094684) were designed for S genes, N genes and complete genome 
amplifications respectively (Table S1). PCR assays were performed with the following 
thermal profile: 95 °C for 5 min, 35 cycles of 95 °C for 30 s, 50 °C for 30 s, and 72 °C for 1 min 
15 s, followed by a final 10 min extension at 72 °C. The products were purified following the 
manufacturer’s instructions of Gel Band Purification Kit (Omega Bio-tek, USA) and then 
cloned into the PMD 19-T vector (Takara, Biotechnology, Dalian, China) and transformed E. 
coli DH5a competent cells. The positive clones were screened out and sent to Beijing 


Genomics Institute (Shenzhen, Guangdong, China) for further sequencing. 


2.4 | Sequence alignment and phylogenetic analysis 


The nucleotide sequences were assembled and aligned using the DNASTAR program 
(DNAStar V7.1, Madison, WI, USA). Phylogenetic trees were constructed using the 
neighbour-joining method in MEGA 7.0 software with bootstrap analysis of 1,000 replicates. 
Percentages of replicate trees in which the associated taxa clustered are shown as nearby 
branches (Chenna et al., 2003; Tamura, Nei & Kumar, 2004; Kumar, Stecher & Tamura, 


2016). 
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3 | RESULTS 


3.1 | SADS-CoV detection 


Of 236 clinical diarrhea samples from 45 swine farms, 53 samples collected from eleven 
swine farms between August 2016 and May 2017 were tested positive for SADS-CoV and 
the locations of these eleven farms were indicated in Figure 1. The farms SM, GL, DL, DE, 
TP, WT, DCD and SC were in Qingyuan, and the farms LS, ZW and YX were in Shaoguan. 
From these eleven SADS-CoV-positive farms, a total of 170 archived diarrheal samples 
including 55 feces and 115 intestinal contents were collected to further investigate diarrheal 
pathogens including SADS-CoV, PEDV, RV, PDCoV, and TGEV. The infection rates of 
these five pathogens in order were 43.5% (74/170), 78.2% (133/170), 21.8% (37/170), 8.8% 
(15/170) and 0% (0/170). SADS-CoV was tested positively both in the intestinal and fecal 
samples and the detection rates were 49.6% (57/115) and 30.9 % (17/55), respectively, which 
were both the second highest among the five tested pathogens. PEDV had the highest positive 
rate at 84.3% (97/115) for intestinal and 65.5% (36/55) for fecal samples. TGEV was tested 


negatively in all 170 diarrheal samples (Table S2). 


The cases of individual infection or co-infection for SADS-CoV, PDCoV, PEDV, and 
RV were also detected. More than half the samples were only infected by one of the four 
viruses at a rate of 58.82% (100/170), and the individual infection of PEDV had the highest 
rate of 37.06% (63/170). In the 74 SADS-CoV-positive samples, 28 samples were infected 


alone by SADS-CoV, and the rest 46 samples were infected by SADS-CoV combined with 
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one to three other pathogens. The co-infection of SADS-CoV and PEDV occurred most 


frequently with the rate of 17.65% (30/170) (Table 1). 


The earliest cases with SADS-CoV-positive infection were found in the farms LS, TP 
and ZW on August 2016 according to Table S2. The intestinal contents from these eleven 
farms were tested positive for SADS-CoV and the positive rates ranged from 25% to 77.78%. 
While the fecal samples from two farms ZW and YX were tested negative for SADS-CoV, 
but positive for PDCoV, PEDV and RV infections. The positive rates of SADS-CoV in the 
fecal samples of other nine farms varied from 13.33% to 100%. The farm LS had the highest 
positive rate of 77.8% (7/9) for intestinal samples and the farms DCD and SC had 100% 


positive rates for fecal samples. 


3.2 | Sequence identities and phylogenetic analyses 


The two complete genomes (accession number MG605090 and MG605091), five N genes 
(accession number MG605087, MG605088, MG605089, MG775251 and MG775253) and 
five S genes (accession number MG605084, MG605085, MG605086, MG775250 and 
MG775252) were sequenced and available in GenBank. The lengths of the complete genome, 
N gene and S gene were 27163bp, 1128bp and 3390bp respectively. The two new complete 
sequences, five N genes and five S genes sequences each shared high sequence identities of 
99.8% - 100%, 99.8% - 100% and 99.6% - 100% with previous reported complete genomes, 
N genes and S genes of SADS-CoV. The two full-length sequences of SADS-CoV showed 
95% - 95.1% nucleotide identities with four bat coronavirus HKU2 strains, and the nucleotide 


identities between the N genes of SADS- CoV studied here and those of bat coronavirus 
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HKU2 strains ranged from 94.2% to 94.3%. The S genes of SADS- CoV in this paper shared 
low nucleotide identities of 79.9% - 80% with the S genes of bat coronavirus HKU2 strains. 
Phylogenetic analyses based on the complete genomes, N genes and S genes indicated that all 
SADS-CoV sequences clustered together to form a well defined branch and group with four 
bat coronavirus HKU2 strains In the phylogenetic trees of the full-length genomes and N 
genes, all SADS-CoV strains were phylogenetically located within the AlphaCoVs, while the 
phylogenetic tree of S genes exhibited that all SADS-CoV sequences, four bat HKU2 and 
two BtRf-AlphaCoV strains separated from other AlpaCoVs sequences and clustered with the 


BeltaCoV sequences (Figures 2, 3, 4). 


4 | DISCUSSION 


In this study, we performed the retrospective investigation of SADS-CoV in 45 pig farms 
from Guangdong Province based on 236 diarrhea samples. Our results showed that the first 
SADS-CoV positive sample was collected in August 2016 from the farm LS with a history of 
diarrhea, as well as from other two farms TP and ZW, which indicates that SADS-CoV has 
emerged in pigs in China at least since August 2016. And this time point is five months 
earlier than the first discovered time reported by our previous study (Zhou et al., 2018). As 
the same time, clinical signs of SADS-CoV during the retrospective investigation included 
sever and acute vomiting and diarrhea, leading to death in piglets that were less than five 


days of age with a mortality rate of around 5096. These clinical presentations were similar to 
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those signs in the large scale outbreak of SADS-CoV reported by Zhou et al. (2018), except 


the mortality rate in piglets later increased to 90%. 


Based on the rates of infection documented in our work, it revealed that PEDV (78.25%) 
was still the primary cause of the porcine diarrhea, which is consistent with previous studies 
that PEDV has been considered to be the major pathogen responsible for the porcine diarrhea 
epidemic in China since 2010 (Sun et al., 2012; Ge et al., 2013; Zhao et al., 20216 ). Except 
PEDV, SADS-CoV had the second high infection rate (43.53%), which exhibits an evidently 
prevailing tendency in pigs. Meanwhile, of the SADS-CoV positive samples, 62.296 were 
co-infections with one to four other viruses, revealing a high prevalence of co-infection in the 
sampled farms. Notably, our results also showed that SADS-CoV and PEDV were 
simultaneously present in all co-infection samples of SADA-CoV. Considering the fact that 
PEDV has caused prior outbreaks of the porcine diarrhea at the pig farm where Zhou et al. 
(2018) later reported the occurrence of SADS-CoV, PEDV seems to be able to contribute to 
the infection of SADS-CoV in pigs. Thus, SADS-CoV may be pathogenic as secondary 
infection following the infection of PEDV, which needs further studies to better understand 


the pathogenesis of this novel coronavirus. 


The phylogenetic relationships of SADS-CoV sequences were also identified in this 
study. The results showed that all SADS-CoV sequences clustered together to form an 
independent branch and separated from other viral sequences in the genus Alphacoronavirus. 
Our results also indicated that both the complete genomes, N genes and S genes of all 
SADS-CoV strains shared the highest nucleotides identifies with those corresponding 
sequences of four bat coronavirus HKU2 strains. In this work, The phylogenetic trees of full 
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length genomes and S genes of SADS-CoV sequences showed that the SADS-CoV branch 
clustered with these four HKU2 strains, which is same to previous results (Gong et al., 2017; 
Pan et al., 2017; Zhou et al., 2018). Besides the genomes and S genes, the tree of N genes in 
our study revealed the identical result too. So far, a total of eight full-length genomes of 
SADS-CoV have been reported in Guangdong Province of China (Gong et al., 2017; Pan et 
al., 2017; Zhou et al., 2018; this study). The two new genomes of SADS-CoV sequences in 
this work shared 100% nucleotides identities with the sequence MF167434 published by 
Gong et al. (2017) and our four previously reported sequences (Zhou et al., 2018), and shared 
99.896 nucleotides identities with the sequence MF370205 studied by Pan et al. (2017). The 
results suggest that these eight SADS-CoV sequences may come from the same origin. Only 
the phylogenetic tree of S genes in our work showed that sequences of the AlphaCoV were 
divided into two sublineages, AlphaCoV1 which contained all SADS-CoV sequences and 
AlphaCoV2, clustering together with sequences of the BeltaCoV and the DelatCov, 
respectively. And this result was consistent with the study of Pan et al. (2017). As a newly 
discovered coronavirus, the availability of SADS-CoV sequences data is limited which 
prevents better understandings of the molecular epidemiology of this virus. Meanwhile, being 
a RNA virus, SADS-CoV may mutate rapidly and exhibit high genetic differences 
(Drummond et al., 2003; Kühnert et al., 2011). So, new more available sequences data are 
warranted to give a deep insight of viral genetic origin, evolution and transmission patterns of 


SADS-CoV. 
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In summary, our retrospective study suggests that SADS-CoV has emerged in pigs in 
China at least since August 2016. The severe clinical symptoms and the strong transmission 
in short term highlights the urgency to develop effective measures to control this new 
discovered virus. Further studies, such as the epidemiology, virology and pathobiology 


should be performed for better understanding of SADS-CoV in China. 
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Figure Legends 


Figure 1. Locations of positive farms for study of SADS-CoV in Guangdong Province, 


China. 


Figure 2. Phylogenetic analysis of the complete genomes of SADS-CoV and reference 
coronavirus species. The tree was constructed using MEGA 7.0 software with 
neighbor-joining methods and 1,000 replicate sets on bootstrap analysis. Two new complete 


genomes studied in this work were indicated with " black solid circles ". 


Figure 3. Phylogenetic analysis of the N genes of SADS-CoV and reference coronavirus 
species. The tree was constructed as per Figure 2 above. Five new sequences of N genes 


studied in this work were indicated with "black solid circles". 
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Figure 4. Phylogenetic analysis of the S genes of SADS-CoV and reference coronavirus 
species. The tree was constructed as per Figure 2 above. Five new sequences of S genes 


studied in this work were indicated with "black solid circles". 


Table 1. Screening information of SADS-CoV, and co-infections of different pathogens 


detected in north Guangdong Province of China between August 2016 and May 2017. 


Table S1. Primers for amplifications for the competed genome, S and N genes of 


SADS-CoV. 


Table S2. Detection of SADS-CoV, PEDV, RV, PDCoV and TGEV in diarrheal samples 


from eleven farms of North Guangdong Province. 
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Table 1 


Four pathogens 








One pathogen positive samples no.(%, Two pathogens positive samples SADS-CoV 
170) (% 170) Three pathogens positive samples no.(%, n=170) positive samples 
n= no.(%, n=17 pus 
Month no.(%, n-170) positive 
s of farms 
year SADS-C PDCo TGE SADS-CoV+ PEDV+R PEDV+PD SADS-CoV+PE SADS-CoV+PEDV SADS-CoV+PEDV+R __ (Positive/T 
PEDV RV PEDV+RV+P tal 
oV V V PEDV V CoV DV+RV +PDCoV V+PDCoV otal) 
DCoV 
Aug, 
3 4 0 0 0 2 3 1 0 0 1 0 
2016 
Sep, 
0 6 0 2 0 0 1 1 0 0 0 0 
2016 
Oct, 
0 2 1 0 0 0 3 0 0 0 0 0 
2016 
Nov, 
1 8 1 1 0 1 4 1 1 1 0 2 11/45 
2016 
Dce, 
2 8 0 1 0 3 0 1 3 0 0 0 
2016 
Jan, 
2 3 2 0 0 10 3 0 2 0 0 0 
2017 
Feb, 
6 7 0 0 0 4 3 1 2 0 0 0 
2017 
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9 0 1 0 4 0 0 4 0 0 0 
2017 
Apr, 
10 10 0 0 0 6 0 0 0 1 0 0 
2017 
May,2 
6 0 0 0 0 1 0 0 0 0 0 
017 
28(16.4 63(37.0 4(23 5(2.9  O(0 18(10.5 
30(17.65%) 5(2.94%) 12(7.06%) 2(1.18%) 1(0.59%) 2(1.18%) 
Total 7%) 6%) 5%) 4%) %) 9%) 





SADS-CoV positive samples : 74(43.53%) 
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63 Bat coronavirus HKU5-1 | NC009020. 1 
id Bat coronavirus HKU9-1 | NC009021.1 NC 009021.1 
Bat SARS coronavirus HKU3-1 | DQ022305.2 
SARS coronavirus | NC004718.3 
56 Bat SARS-like coronavirus RsSHC014 | KC881005.1 
55 L Bat SARS-like coronavirus WIV1 | KF367457.1 
53 Avian infectious bronchitis | NC001451 
European turkey coronavirus 080385d | KR822424.1 


70 


67 





Pa 








68 Bulbul coronavirus HKU 11-796 | FJ376620.1 
52 [ Porcine deltacoronavirus isolate PDCoV/CHJXNI2/2015 | KR131621.1 
51 | Porcine coronavirus HKU15 strain HKU15-44 | NC016990.2 
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AplhaCoV 


BetaCoV 


GammaCoV 


DetalCoV 


Porcine enteric alphacoronavirus GDS04 | MF167434.1 
@ SADS-CoV/CN/GDTP/2017 N gene | MG605088 
@ SADS-CoV/CN/GDZW/2017 N gene | MG605089 
@ SADS-CoV/CN/GDLS/2017 N gene | MG775253 
@ SADS-CoV/CN/GDGL/2017 N gene | MG775251 
@ SADS-CoV/CN/GDDL/2017 N gene | MG605087 
Rhinolophus bat coronavirus HKU2 CH/GD-01/2017/P2 | MF370205.1 
SADS-CoV_FarmD | MF094684.1 
100| | SADS-CoV_FarmC | MF094683.1 
SADS-CoV_FarmB | MF094682.1 
SADS-CoV_FarmA | MF094681.1 
Bat coronavirus HKU2 | NC009988.1 
100 Bat coronavirus HKU2 strain HKU2/HK/298/2006 | EF203066.1 
H Bat coronavirus HKU2 strain HKU2/HK/33/2006 | EF203067.1 
Bat coronavirus HKU2 strain HKU2/HK/46/2006 | EF203065.1 
BtRf-AlphaCoV/YN2012 | KJ473808.1 
100 ' BtRf-AlphaCoV/YN2012 | NC028824.1 
41 Bat coronavirus 1A | NC010437.1 
94 Bat coronavirus HKU8 strain AFCD77 | EU420139.1 
Rousettus bat coronavirus HKU10 | NC01887.1 
Bat coronavirus CDPHE15/USA/2006 | NC022103.1 
96 Scotophilus bat coronavirus 512 | N009657.1 


67 E Porcine epidemic diarrhea virus strain GDSO1 | KM089829.1 


100 





29 
AplhaCoV 


100 |r Porcine epidemic diarrhea virus strain CV777 | KT323979.1 
72L Porcine epidemic diarrhea virus| NC003436.1 
99 Human Coronavirus NL63 | NC005831.2 
100 NL63-related bat coronavirus strain BtKYNL63-9a | NC032107.1 
38 229E-related bat coronavirus strain BEKY229E-1 | KYO73747.1 
Camel alphacoronavirus isolate Riyadh/Ry141/2015 | NC028752.1 
100 Human coronavirus 229E | NC002645.1 
TGEV Purdue P115 | DQ811788.1 
100 | | PRCV ISU-1 | DQ811787.1 
100' Transmissible gastroenteritis virus | AJ271965.2 


100 


91 
54 - Bovine coronavirus | NC003045.1 


83 Porcine hemagglutinating encephalomyelitis virus | NC007732.1 
100| | PHEV CoV USA-15TOSU1582 | KY419113.1 
100 Human coronavirus OC43 strain SC2481 | KY983583.1 
Mouse hepatitis virus strain MHV-A59 C12 mutant | NC001846.1 
100 Murine hepatitis virus strain JHM complete genome | AC000192.1 
Bat coronavirus HKU9-1 | NC009021.1 NC009021.1 
100 Coronavirus Neoromicia/PML-PHE 1/RSA/2011 | KC869678.4 
38 Middle East respiratory syndrome coronavirus | NC019843.3 
B 100 Bat coronavirus HKU5-1 | NC009020.1 
Bat coronavirus HKU4-1 | NC009019.1 
55 Bat SARS coronavirus HKU3-1 | DQ022305.2 
100| p Bat SARS-like coronavirus RSSHCO14 | KC881005.1 
99 | Bat SARS-like coronavirus WIV1 | KF367457.1 
38- SARS SARS ocoronavirus | NC004718.3 
European turkey coronavirus 080385d | KR822424.1.1 
Bulbul coronavirus HKU11-796 | FJ376620.1 
Porcine deltacoronavirus isolate PDCoV/CHJXNI2/2015 | KR131621.1 DetalCoV 
99 L Porcine coronavirus HKU15 strain HKU15-44 | NC016990.2 


96 
BetaCoV 


| GammaCoV 


100 


0.1 
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Rhinolophus bat coronavirus HKU2 CH/GD-01/2017/P2 | MF370205.1 
6q @ SADS-CoV/CN/GDTP/2017 S gene | MG605085 


300 @ SADS-CoV/CN/GDGL/2017 S gene | MG775250 


Porcine enteric alphacoronavirus GDS04 | MF167434.1 
@ SADS-CoV/CN/GDZW/2017 S gene | MG605086 
@ SADS-CoV/CN/GDLS/2017 S gene | MG775252 
97| @ SADS-CoV/CN/GDDL/2017 S gene | MG605084 
96 SADS-CoV FarmC | MF094683.1 
SADS-CoV FarmB | MF094682.1 
SADS-CoV FarmA | MF094681.1 
SADS-CoV FarmD | MF094684.1 
Bat coronavirus HKU2 | NC009988.1 
Bat coronavirus HKU2 strain HKU2/HK/46/2006 | EF203065.1 
00 | Bat coronavirus HKU2 strain HKU2/HK/33/2006 | EF203067.1 
Bat coronavirus HKU2 strain HKU2/HK/298/2006 | EF203066.1 
BtRf-AlphaCoV/YN2012 | KJ473808.1 
100 ! BtRf-AlphaCoV/YN2012 | NC028824.1 
Bat coronavirus HKU9-1 | NC009021.1 
100 [— Bovine coronavirus | NC003045. 1 
100 Human coronavirus OC43 strain SC2481 | KY983583.1 
PHEV CoV USA-15TOSU1582 | KY419113.1 
90 100 100 — Porcine hemagglutinating encephalomyelitis virus | NC007732.1 
Mouse hepatitis virus strain MHV-A59 C12 mutant | NC001846.1 
100 Murine hepatitis virus strain JHM complete genome | AC000192.1 
99 Bat coronavirus HKU4-1 | NC009019.1 
100 Bat coronavirus HKU5-1 | NC009020.1 
Coronavirus Neoromicia/PML-PHE1/RSA/2011 | KC869678.4 
95 Middle East respiratory syndrome coronavirus | NC019843.3 
d Bat SARS coronavirus HKU3-1 | DQ022305.2 
SARS coronavirus | NC004718.3 
100 Bat SARS-like coronavirus RSSHC014 | KC881005.1 
100 t Bat SARS-like coronavirus WIV1 | KF367457.1 
European turkey coronavirus 080385d | KR822424.1.1 | GammaCoV 
100 r Porcine deltacoronavirus isolate PDCoV/CHJXNI2/2015 | KR131621.1 
100 Porcine coronavirus HKU15 strain HKU15-44 | NC016990.2 DetalCoV 
Bulbul coronavirus HKU11-796 | FJ376620.1 
64 Bat coronavirus HKU8 strain AFCD77 | EU420139.1 
85 Rousettus bat coronavirus HKU10 | NC01887.1 
100 Bat coronavirus 1A | NC010437.1 
100 [— Camel alphacoronavirus isolate Riyadh/Ry141/2015 | NC028752.1 
100 Human coronavirus 229E | NC002645.1 
100 229E-related bat coronavirus strain BUKY229E-1 | KY073747.1 
3 Human Coronavirus NL63 | NC005831.2 
NL63-related bat coronavirus strain BtKYNL63-9a | NC032107.1 AplhaCoV2 







AplhaCoV1 


100 





95 


BetaCoV 
61 





100 





98 


100 | Transmissible gastroenteritis virus | AJ271965.2 
31 100 | ' TGEV Purdue P115 | DQ811788.1 
L PRCV ISU-1 | DQ811787.1 
90 Bat coronavirus CDPHE15/USA/2006 | NC022103.1 
Scotophilus bat coronavirus 512 | NC009657.1 
Porcine epidemic diarrhea virus strain GDSO1 | KM089829.1 
Porcine epidemic diarrhea virus strain CV777 | KT323979.1 





83 


100 
97 L Porcine epidemic diarrhea virus] NC003436.1 
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